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Abstract 

In this paper, we present a three-mode dynamo model which describes both supercritical and 
subcritical dynamo transitions. The nature of dynamo transition changes from supercritical to 
subcritical as the magnetic Prandtl number is decreased, consistent with the numerical results 
in the spherical-shell and the Taylor-Green dynamo. We also perform a detailed analysis of the 
hysteresis zone of a subcritical dynamo using direct numerical simulations. Numerical simulation 
and model analysis show that the sets of initial conditions, called the basin of attraction, of the 
no-dynamo and the dynamo states are separated by an unstable manifold. 
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I. INTRODUCTION 



Many natural and engineering systems exhibit transitions from one state to another. Some 
of the prominent examples of transitions are water to vapor, paramagnetic to ferromagnetic, 
conduction to convection (e.g, in Rayleigh Benard convection), laminar to turbulent flow 
in channels, etc. Some of these transitions are continuous, i.e., the order parameter grows 
smoothly from zero as the control parameter is increased, while some others exhibit a dis- 
continuity or a finite jump in the order parameter. The former class of transitions is called 
supercritical, while the latter is called subcritical. 

Magnetic field generation or dynamo process in astrophysical systems exhibits transition 
from no-dynamo or fluid state to a magnetic state l|, This transformation is referred 
to as "dynamo transition". Laboratory experiments report supercritical transition js-G], 
but numerical simulation and models exhibit both supercritical |7H9| and subcritical transi- 
tions 

lel-Q- In this paper, we explain these transitions using a low- dimensional model and 
relate them to the strength of the nonlinear interactions. 

Nature of dynamo transition can play a very crucial role in our understanding of astro 
and planetary dynamos. For instance. Mars is believed to have no working dynamo, but it 
exhibits strong crustal magnetism. Kuang et al. 15| conjectured that the sudden termination 
of the Martian dynamo could be due to a subcritical dynamo transition. Christensen et 
al. 0] and Morin and Dormy 16] used numerical simulations to study dynamo mechanism in 
rotating spherical-shells; they observed both supercritical and subcritical dynamo transitions 
in such systems. The control parameters for dynamo simulations are the magnetic Prandtl 
number Pm = v/rj, the Reynolds number Re = UL/v, and the magnetic Reynolds number 
Rm = UL/rj, where v^rj are the kinematic viscosity and magnetic diffusivity respectively, 
and f/, L are the large-scale velocity and length scales. These studies also indicate that 
the dynamo transition changes from supercritical to subcritical as the magnetic Prandtl 
number Pm is decreased. Similar dependence on Pm has been observed in the Taylor-Green 
dynamo j^, [lo|, [l^ and shell models [l^. Based on geometrical similarities and the 
aforementioned Pm dependence, Yadav et al. [l^ conjectured a connection between the 
Taylor-Green dynamo and the spherical-shell dynamo. 

Low-order dynamo models, despite of their simplicity, can nonetheless provide very in- 
teresting insights. BuUard 17|] used a homopolar disk dynamo to study the basic dynamo 
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mechanism, and, subsequently, Rikitake [l8| used a coupled disk dynamo model to quali- 
tatively discern the geomag net. ..eve.a.s P*as and Fauv. Q constructed a >o„-oMc. 
dynamo model based on amplitude equations involving dipolar and qudrupolar modes to 
study magnetic field reversals in the Von-Karman-Sodium (VKS) experiment s], [2o|. Later 
Gissinger et al. {21 1 proposed a three- mode model for similar analysis. Verma et al. 22 1 
constructed a six-mode model to study the importance of helicity in dynamo transitions; 
they report a supercritical pitchfork dynamo transition. A subcritical dynamo transition has 
been demonstrated by Fedotov et al. 23|] in a perturbative af2-dynamo model. Krstulovic 



et al. 



9| devised a model containing two nonlinear equations to explain the Pm dependence 



of the supercritical and subcritical dynamo transition. Weiss [2J] provides a detailed review 
of the low-order dynamo models which are related to the planetary and solar dynamos. 

In this paper we construct a three-mode model that captures both supercritical and 
subcritical dynamo transition. The model shows a very good agreement with the recently 
performed numerical simulations of the Taylor- Green dynamo jsi [l^. 

The paper is structured in the following manner. In Sec. [Tll we discuss one-dimensional 
equations which show supercritical and subcritical bifurcations. Section UTTl contains a simu- 
lation study of the Taylor-Green dynamo in a subcritical regime. We propose a three-mode 
model in Sec. HVl We conclude in Sec. |V] 



II. ONE-DIMENSIONAL MODEL 



Our three-mode dynamo model, to be described later, is motivated by a one- dimensional 



(ID) model that exhibits supercritical and subcritical transitions 
the model 
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26j . The first part of 



X = CiX + C-sX^ 



(1) 



exhibits supercritical transition. For negative values of C3, the solution X = is stable for 
Ci < 0. For Ci > 0, this solution becomes unstable, and a new stable solution ■>/— C1/C3 
is born. The bifurcation diagram for this system with C3 = — 1 is illustrated in Fig. [TJ^a). 
Note that the new solution exists only when C3 < 0; for Ci,C3 > 0, the above system is 
unstable. 

When C3 > 0, a higher-order term C^X^ with C5 < is added to stabilize the system. 
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i.e., 



X = CiX + CsX^ + C^X\ (2) 

Fig. [U^b) illustrates the bifurcation diagram of this system for C3 = 1 and C5 = — 1. It 
exhibits a subcritical pitchfork bifurcation and a jump in the value of X at the bifurcation 
point Ci = 0. In addition, the system undergoes a saddle-node bifurcation at Ci = —0.25. 
The hysteresis cycle represented by the arrows in Fig. [H^b) results due to these bifurcations. 
The hysteresis zone has multiple final states, as evident from Fig. [T]^b). For —0.25 < Ci < 0, 
initial conditions below the unstable manifold (solid curve connecting hollow red-dots) settle 
down to X = states (filled blue-dots on X = axis), while the ones above the unstable 
manifold settle down to X > states (filled blue-dots). A similar situation exists for the 
X < states. Thus, initial condition plays a critical role in determining the final state of 
a subcritical system. In the following discussion, we will construct a three-mode dynamo 
model that exhibits similar features as the aforementioned ID model. 

In the next section we report numerical results of a Taylor-Green dynamo in a subcritical 
regime. We show that the numerical results are in good agreement with the behavior of ID 
model described by Eq. (|2]). 



III. SUBCRITICAL DYNAMO TRANSITION IN TAYLOR-GREEN DYNAMO 

The fluid velocity u and the magnetic field B in dynamo mechanism are governed by the 
magnetohydrodynamic (MHD) equations: 

dtu + (u ■ V)u = -Vp + (J X B) + uV^u + F, (3) 

dtB = V X (u X B) + r/V^B, (4) 

V ■ u = 0, (5) 

V ■ B = 0, (6) 

where J, F, p, u, and rj represent current density, external forcing, hydrodynamic pressure, 
kinematic viscosity, and magnetic diffusivity respectively. Note that J = V x B. The density 
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FIG. 1: (Color online) Bifurcation diagrams for Eq. ([T|) and Eq. ([2]) portraying (a) supercritical 
and (b) subcritical transitions. The filled-blue and the hollow-red dots represent the stable and 
unstable states. In figure (b), the arrows depict the hysteresis cycle of the subcritical transition. 



of the fluid is chosen to be unity. The Taylor-Green (TG) forcing, 



sm 



.{kox) cos{kQy) cos{koz) 
— cos(A;oa^) sin(A;oy) cos^koz) 







(7) 



is one of the popular forcing schemes to study dynamo transition jsijiol. [l^. Here, Fq is the 
forcing amplitude, and ko (=2) defines the length scale of the forcing. Yadav et al. 
showed that dynamos with the TG forcing exhibit supercritical and subcritical dynamo 



transition for Pm = 1 and Pm = 0.5 respective^ 



the bifurcation reported by Yadav et al. 



lower than 0.5 could also be subcritical [id, 12|. Near the dynamo transition, the nature of 



y. It has been shown that dynamos with Pm 



14l | is similar to the aforementioned ID model. 




FIG. 2: (Color online) Bifurcation diagram constructed using various dynamo states of the TG 
dynamo. The filled blue-stars connected by solid line represent a stable dynamo branch. The 
runs with hollow blue-stars as initial conditions go back to the dynamo branch (filled blue-stars), 
while the runs with hollow red-triangles as initial conditions settle down to the no-dynamo states 
(filled red-triangles). The initial conditions denoted by hollow purple-squares near the hysteresis 
boundary stabilize to a different dynamo branch represented by filled purple-squares. 



The reader is referred to the Fig. 4 of Ref. 
Pm = 0.5), which are similar to Fig. [1] 



(for Pm = 1) and the Fig. 3 of Ref. [ij (for 



Subcriticality in the dynamo transition has been inferred based on a jump in the mag- 
netic energy and a hysteresis near the transition jlo|, [l^ |2^. The generic subcritical bifur- 
cation described by Eq. ([2]) exhibits more features, e.g., existence of subcritical-pitchfork 
and saddle-node bifurcations. Also, in the hysteresis window, the system stabilizes to either 
X = or X 7^ depending on the initial condition. Using direct numerical simulations 
(DNS) we explore whether these rich features are present in the subcritical regime of the 



TG dynamo or not. We use a pseudospectral code Tarang [28| to numerically solve Eqs. ([3])- 
with Pm = 0.5 {u = 0.1 and t] = 0.2) in a cube of dimensions {2tt)^ having periodic 
boundary conditions on all sides. We discretize the simulation box using 64'^ grid points. 
Interested reader is referred to Ref. [l^ for more details about the numerical procedure. 
Figure [2] contains results obtained using DNS of the TG dynamo for Pm = 0.5. As 
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reported by Yadav et al. [l^, at Fq = 15.8 we first observe a fixed-point dynamo state 
(temporally non-fiuctuating); here, the magnetic energy shows a finite jump to ~ 2.6 
(the topmost star of the filled blue-star trail). However, when the dynamo state at Fq = 15.8 
is used as an initial condition and the forcing amplitude is gradually decreased, the dynamo 
state continues along the trail till Fq ^ 15.15, at which point there is a sudden dip to the 
fiuid state (filled red-triangles). Thus, Fq = 15.15 : 15.8 is the hysteresis band for Pm = 0.5 
TG dynamo. To explore the basins of attraction of the E^ = and the E^ > states, we 
construct a set of initial conditions using the dynamo states on the filled blue-star trail. 

For the initial conditions, we keep the velocity field same as that of the dynamo state, 
but the magnetic field is quenched by a factor, e.g., 2, 4, 6, etc. The set of initial conditions 
used in this work are indicated by hollow data points in Fig. [21 Interestingly, for a given 
Fo, all the runs with hollow blue-stars as initial conditions stabilize to the corresponding 
dynamo states (filled blue-stars). On the contrary, the runs with hollow red-triangles as 
initial conditions settle to the corresponding no-dynamo states (filled red-triangles). Thus, 
in the hysteresis zone, the regions with hollow blue-stars and hollow red-triangles form the 
basins of attraction of the dynamo and no-dynamo states respectively. The dashed black 
curve separates the two regions, a feature very similar to that observed for the subcritical 
system of Eq. ([2]). 

An another feature appears in our simulation: the runs with the initial conditions denoted 
by hollow purple-squares stabilize to the corresponding filled purple-square dynamo state 
(at Fq = 15.15), which belongs to an another dynamo branch. Note that this new branch 
originates from the edge of the hysteresis, and it is related to the high- dimensionality of 
the system. The other dynamo states on the purple trail were constructed using the filled 
purple-square dynamo state for Fq = 15.15. The above analysis demonstrates that the 
system is very sensitive to the initial condition, and we need to keep this in mind while 
studying dynamo. 

The TG dynamo under investigation has many degrees of freedom. Yet, similarities be- 
tween the bifurcation diagrams of Fig. lU^b) and Fig. [2] indicate that the system dynamics 
could be understood using the important modes of the system. Yadav et al. ^, |l4| showed 
that for the TG forcing with ko = 2, the most dominant velocity and magnetic Fourier modes 
are u(2, 2, 2) and B(0, 0, 1) respectively, and they couple via B(— 2, —2, —1) Fourier mode. 
The dominance of the large-scale modes is consistent with the results of the VKS experi- 



7 



ment 



201 ] and the planetary dynamo simulations |7|] in which the dipolar and quadrupolar 



magnetic field modes dominate near the transition. Motivated by these observations, we 
construct a three-mode model. We will show in the next section that this three-mode model 
captures the transition from supercritical to subcritical dynamo quite nicely. 

IV. THREE-MODE DYNAMO MODEL 



In our model, we consider a velocity mode u and two magnetic modes Bi and B2 which 
are coupled nonlinearly with each other. The equations of the model are 



ii = f — vu — B1B2, 
Bi = auB2 — rjBi, 
B2 = uBi - 7]B2 + ^uBf, 



(8) 
(9) 
(10) 



where u and rj are the kinematic viscosity and magnetic diffusivity respectively, / is the 
external force acting on the velocity field, and a, /3 are positive constants. The nature of 
the nonlinear interactions resembles the MHD interactions where dtB{]<.) ~ u{p)B{q) and 
dtu{k.) ~ B{p)B{q) with k = p + q. Also, our model respects the symmetry of the MHD 
equations under B — t- — B transformation. The coefficients of the nonlinear terms UB2 and 
uBi are different since Bi and B2 correspond to different wavenumbers. The cubic term 
uBf, induced by a higher order of nonlinearity, plays an important role in dynamo transition 



2l| except 



of this model. Our model has certain similarities with that of Gissinger et al. 
the cubic interaction. 

For the aforementioned model, the first dynamo state is a fixed point, which we study 
perturbatively. We compute fixed point solutions by setting the right-hand-side of Eqs. ([8]) 
and ( ITOj) to zero, which yields 

uBi + l3uBl 



Bo 



V 



f 

u = — 

V 



1 + 



B\(\^fiBX, 

VT] 



-1 



(11) 



(12) 



Substitution of the above expressions of B2 and u in Eq. ([9]) provides 



5i 



{Bi + /3B!) 



B!{l + f3Bl 

VT] 



vBi 



(13) 



An approximation of Eq. (fT3l) to the fifth order in Bi yields 



Bi 



Bi 



^ HP - 2) 



Bl 



ap 3 



-4/3) 



Bl 



(14) 



Ci Ca 

Note that Eq. f lT^ has the same form as Eq. ([2]), and the nature of its transition is governed 
by the coefficients C3 and C5. 

The dissipative parameters v and 77 take positive values. In addition, we assume that 
Q!,/3 > 0. Under these assumptions, the following properties can be deduced for Ci,C3, and 
C5: 

1. For a given values of v and r/, Ci changes sign depending on the force parameter /, 
hence / is the control parameter. 

2. The sign of C3 depends on the term {vrjP — 2). For a given z/, C3 < for rj < rjc = 
2/(/3i/), which provides a supercritical transition. For rj > t]c or uriP > 2, we observe 
subcritical transition. We can deduce the critical magnetic Prandtl number Pmc that 
separates the supercritical and subcritical dynamo using the condition (3ut]c = 2, which 
yields 



Pmc = — 

Vc 



(15) 



3. The sign of C5 is determined by (| — 1^77 /3). For subcriticality, we require C3 > and 
C5 < (see Sec. [Tll), which are satisfied by the condition urjP > 2. 

Thus the structure of Eq. (fT^ is same as that of one-dimensional model discussed in Sec. [Tll 
We construct a potential function for Bi using Eq. ( 11^ as 



V{B^ 



BidBi 



2 ^ A ^6 



(16) 
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26|. 



This function provides an elegant way of visualizing the nature of the dynamo states 
In this description, the stable and the unstable dynamo states are represented by "valleys" 
and "hills" , respectively, of the potential function. 

In the following discussion we will illustrate the supercritical and subcritical dynamo 
transitions exhibited by the three-mode model. 
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A. Supercritical Transition 



Here we present a representative example of a supercritical dynamo by taking v = 10, 
?7 = 1.5, /3 = 0.1, a = 2 for which C3 < 0. For these parameters, Pm^ = 5 and Pm = 
10/1.5 ~ 6.6. Here, dynamo transition takes place at 



Figure [3](a) presents the potential functions, which exhibit a valley at Si = for / < /c, 
and two valleys at i?i 7^ (labeled by blue dots) for / > fc. Hence, the no-dynamo states 
are stable for / < f^, and the dynamo states are stable for f > fc- A bifurcation diagram 
constructed using these states is illustrated in Fig. [3](b); filled-blue (hollow-red) dots are the 
stable (unstable) states. 

B. Subcritical Transition 

To illustrate a subcritical dynamo transition in the three-mode model, we take 77 = 3, 
along with i/ = 10, /3 = 0.1, a = 2. Here Pm = 10/3 ~ 3.3. These parameters yield C3 > 
and C5 < 0, thus yielding a subcritical dynamo transition. The critical forcing fc for this 
set of parameters is 30/-\/2 ^ 21.2. Near / = fc, we plot the potential functions and the 
bifurcation diagram in Fig. Hl^a) and (b) respectively. At f = fc we observe a finite jump in 
Bi (see Fig. ID^b)). As we decrease / from fc with the above initial condition, the system 
continues to have nonzero Bi until / = /c ~ 20.94, after which it attains -Bi = state. 

In the potential diagram of Fig. SJ^a), we observe a single valley (stable state) at -Bi = 
for / < and two valleys at ±i?i and a hill at -Bi = for / > fc- For < / < fc, the 
system has three valleys and two hills. Depending on the initial condition, the system can 
"roll" to one of the valleys which explains the sensitivity to initial condition. 

To better understand the role of initial condition on the final state of the system, we plot 
a separate potential function for / = 21.04 in Fig. [5l This forcing value lies in between 
and fc- If the initial condition lies on the dashed region of the potential function, then the 
system will roll down to the Bi = valley (no-dynamo state). On the other hand, if the 
initial condition lies on the solid regions of the potential function, then the system will roll 
to either the Bi > valley or the -Bi < valley, as indicated by the arrows; in these cases. 



/. = ^ - 10.6. 




(17) 
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FIG. 3: (Color online) Plots of (a) Potential function, and (b) the corresponding supercritical 
transition of the dynamical system described by Eq. (I14p . We set i^, ?], a, and /3 equal to 10, 1.5, 
2, and 0.1 respectively. In the top panel, the forcing / is varied from 10.55 to 10.65 in steps of 
0.01. Filled-blue (hollow-red) dots represent the dynamically stable (unstable) states. 



the system will be a dynamo state. This behavior is applicable only for < / < /c- For 
/ > /c, the two Bi ^ unstable hills merge with Bi = state. 

The above results indicate that the three-mode dynamo model captures the supercritical 
to subcritical dynamo transition nicely. Moreover, this change in the nature of the dynamo 
transition is exhibited when we decrease the magnetic Prandtl number. 
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FIG. 4: (Color online) Plots of (a) Potential function, and (b) the corresponding subcritical tran- 
sition described by Eq. (I14p . We set v, r], a, and /3 equal to 10, 3, 2, and 0.1 respectively. In the 
top panel, the forcing / is varied from 20.84 to 21.2 in steps of 0.04. The filled-blue (hollow-red) 
dots represent the dynamically stable (unstable) states in both the figures. 



V. CONCLUSIONS AND DISCUSSIONS 



In this paper, we construct a three-mode dynamo model motivated by the symmetries 
and the nature of the nonlinearities of the MHD equations. It contains a crucial third-order 
term Bf, a representative of a higher-order nonlinearity of the system. Our model exhibits 
both supercritical and subcritical dynamo transitions, with an interesting changeover from 
supercriticality to subcriticality as the magnetic Prandtl number is decreased. The above 
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FIG. 5: (Color online) A plot of the potential function for / = 21.04. Here /c < / < fc- The dashed 
curve represents the initial conditions for the no-dynamo state (filled blue-dot at Bi =0), while 
the solid curve represents the initial conditions for the dynamo state (filled blue-dots at Bi ^ 0). 
The arrows represent the rolling direction of the dynamical system. 



features are consistent with the numerical results of the Taylor-Green 



shell dynamos 



m 



and spherical- 



16| . VKS experiment, which exhibits supercritical dynamo for a very small 
Pm, appears to be contrary to the above behavior. The reason for this disagreement is 
not apparent at present; it could possibly be due to some stabilizing feature of cylindrical 
geometry employed in the experiment. 

In our three-mode model, the dynamo transition depends crucially on the extent of stabi- 
lizing and destabilizing terms. Supercritical to subcritical transition takes place when the Bf 
term of Eq. (fT^ destabilizes the magnetic field, which is subsequently stabilized by an even 
higher order interaction of the equation. This feature is similar to the recent findings of 

n 

Sreenivasan and Jones [27| for subcritical spherical-shell dynamo in which the Lorentz force, 
which otherwise saturates the magnetic field, destabilizes it. 

We also present a detailed analysis of the hysteresis region of a subcritical dynamo tran- 
sition using direct numerical simulation and three-mode model. We show that the final 
state of a subcritical dynamo depends quite critically on the initial condition. The basin of 
attraction of the dynamo and the no-dynamo states are separated by a well-defined unstable 
manifold. 
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The potential function described in this paper has some similarities with those observed 
in the Maxwell's construction in equilibrium thermodynamics as well as phenomenon like 
supercooling. It may be interesting if some features like initial condition sensitivity and 
multiple states discussed in our paper could be observed in such systems for some exceptional 
cases. 

In summary, our dynamo model exhibits supercritical and subcritical dynamo transitions. 
The model also shows a changeover from supercriticality to subcriticality when the magnetic 
Prandtl number is decreased. It would be interesting to make quantitative comparisons of 
the model with realistic dynamos. However, the analysis of the large number of modes in 
such systems is very complex, which may be possible in future studies. 
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